One and two-center processes in high-order harmonic generation in diatomic 
molecules: influence of the internuclear separation 
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We analyze the influence of different recombination scenarios, involving one or two centers, on 
high-order harmonic generation (HHG) in diatomic molecules, for different values of the internuclear 
separation. We work within the strong-field approximation, and employ modified saddle-point 
equations, in which the structure of the molecule is incorporated. We find that the two-center 
interference patterns, attributed to high-order harmonic emission at spatially separated centers, are 
formed by the quantum interference of the orbits starting at a center Cj and finishing at a different 
center C„ in the molecule with those starting and ending at a same center Cj. Within our framework, 
we also show that contributions starting at different centers exhibit different orders of magnitude, 
due to the influence of additional potential-energy shifts. This holds even for small internuclear 
distances. Similar results can also be obtained by considering single-atom saddle-point equations 
and an adequate choice of molecular prefactors. 
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I. INTRODUCTION 

Molecules in strong laser fields have attracted a great 
deal of attention in the past few years. Indeed, strong- 
field phenomena, such as high-order harmonic genera- 
tion, above-threshold ionization, and nonsequential dou- 
ble ionization, may be used as tools for measuring and 
even controlling dynamic processes in such systems with 
attosecond precision This is a direct consequence of 
the fact that the physical mechanisms behind such phe- 
nomena take place within a fraction of the period of the 
laser field. For a typical, titanium-sapphire laser used in 
experiments, whose period is of the order r ~ 2.7fs, this 
means hundreds of attoseconds. 

Explicitly, such phenomena can be described as the 
laser-assisted rescattering or recombination of an elec- 
tron with its parent ion, or molecule 1. At an instant 
t', this electron reaches the continuum through tunnel- 
ing or multiphoton ionization. Subsequently, it propa- 
gates in the continuum, being accelerated by the exter- 
nal field. Finally, it is driven back towards its parent ion, 
or molecule, with which it rccombines or rescatters at a 
later instant t. In the former case, the electron kinetic 
energy is converted in a high-energy, XUV photon, and 
high-order harmonic generation (HHG) takes place [3|. 
In the latter case, one may distinguish two specific sce- 
narios: The electron may suffer an elastic collision, which 
will lead to high-order above-threshold ionization (ATI) 
IH, or transfer part of its kinetic energy to the core, and 
release other electrons. Hence, laser-induced nonsequen- 
tial double (NSDI), or multiple ionization (NSMI) will 
occur. 

For molecules, there exist at least two centers with 
which the electron may recombine or rescattcr. This 
leads to interference patterns which are due to photo- 
electron or high-harmonic emission at spatially separated 
centers, and which contain information about its specific 
structure. In the simplest case of diatomic molecules, 
such patterns have been described as the microscopic 




counterpart of a double-slit experiment 0, Q . 

A legitimate question is what sets of electron orbits are 
most relevant for the two- or many-center interference 
patterns. To understand this issue is a first step towards 
controlling such processes by, for instance, an adequate 
choice of the shape and polarization of the external field. 
In the specific case of diatomic molecules, the electron 
may start and return to the same center Cj , or leave from 
a center Cj and return to a center C v , v ^ j(j = 1,2). 
Hence, in total, there exist four possible processes that 
contribute to the yield. Recently, these processes have 
been addressed in several studies, for above-threshold 
ionization 
ation m, LL 



1 01 . Illl ]. high-order harmonic gener- 
131 . 1141 ] and nonsequential double ionization 
|15|. The vast majority of these studies has been per- 
formed using semi-analytical methods, in the context of 
the strong-field approximation. In this framework, the 
transition amplitude can be written as a multiple inte- 
gral with a slowly varying prefactor and a semiclassical 
action. The structure of the molecule may be either in- 
corporated in the former (Sjlljll [H HJ HH, or in 
the latter @, fl GH, G3, uWuj[ (U On a more spe- 
cific level, when solving these integrals employing saddle- 
point methods, it is possible to draw a space-time picture 
of the laser-assisted rescattering or recombination process 
in question, and establish a direct connection to the or- 
bits of a classical electron in a strong laser field [22| . By 
incorporating the structure of the molecule in the action, 
one obtains modified saddle-point equations which gives 
the one-or two-center scenarios. 

In a previous publication [l4[ , we have addressed this 
issue to a large extent for high-order harmonic genera- 
tion, within the Strong-Field Approximation (SFA). Our 
results suggested that the maxima and minima observed 
in the spectra were due to the quantum interference of 
the processes in which the electron leaves and returns to 
a specific center Cj in the molecule with those in which 
it leaves from Cj, but returns to a different center C„. 
There exist, however, a few ambiguities as far as the in- 
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terpretation of our findings is concerned. For instance, in 
the length-gauge formulation of the SFA, we found addi- 
tional potential energy shifts, which depend on the field 
strength E(t) and in the internuclear separation R. These 
shifts led to a strong suppression of tunnel ionization at 
one of the centers. This could have led to the conclusion 
that the interference between other processes were not 
relevant for the patterns in the spectra. 

In this proceeding, we investigate the role of the one 
and two-center recombination scenarios in more detail. 
In particular, we analyze the above-mentioned poten- 
tial energy shifts and their influence on the spectra, 
for smaller internuclear distances than those taken in 
We also provide an alternative interpretation of 
the results encountered, based on effective prefactors and 
single-atom saddle-point equations. 

This paper is organized as follows. In Sec. UH we 
briefly recall the strong-field approximation HHG transi- 
tion amplitudes. Thereby, we consider the situation for 
which the structure of the molecule is either incorporated 
in the prefactor (Sec. Ill BH . or in the semiclassical action 
(Sec. Ill C[) . Subsequently (Sec. IIII[) , we analyze the role 
of the different scenarios, involving one and two centers, 
in the high-harmonic spectra, either solving the modi- 
fied saddle-point equations (Sec. IIII Aft , or mimicking 
the quantum-interference between different sets of orbits 
by an adequate choice of prefactors (Sec. IIII Bp . Finally, 
in Sec. IIVI we outline the main conclusions of this work. 



II. TRANSITION AMPLITUDES 



with the action 



s(t,t',n,k) = -±J\k+A(T)] 2 dT-i p 



(t-t')+nt (4) 



k(f ) fl^ (*') |^o) ■ Thereby r, e 2 



and the prefactors d rec (k(i)) = (k(t) r.e x |?/>o) and 

d ion {k(t')) = 

H- m t(t'), I p , and O give the dipole operator, the laser- 
polarization vector, the interaction with the field, the 
ionization potential, and the harmonic frequency, respec- 
tively. The explicit expressions for k(t) are gauge depen- 
dent, and will be provided below. Physically, Eq. (|3|) 
describes a process in which an electron, initially in a 
field- free bound-state \ipo), is coupled to a Volkov state 

k(<')^) by the interaction H mt (t') of the system with the 

field. Thereafter, it propagates in the continuum and is 
driven back towards its parent ion, or molecule. At a 
time t, it recombines, emitting high-harmonic radiation 
of frequency fl. 

The above-stated transition amplitude may be ei- 
ther solved numerically, or employing saddle-point equa- 
tions. In this work, we employ the latter method 
and the specific uniform approximation discussed in 
Ref. [25l |. Explicitly, these equations are given by 
the condition that the semiclassical action be station- 
ary, i.e., that d t S(t,t',Q,k) = d v S(t,t',fl,k) = and 
dkS{t,t',n,k) = 0. 

For a single atom placed at the origin of the coordinate 
system, this leads to 



A. General expressions 



[k + A(t')] =-2I p 



(5) 



As a starting point, we will underline our main as- 
sumptions with regard to the diatomic bound-state wave 
functions. We consider frozen nuclei, the linear combi- 
nation of atomic orbitals (LCAO) approximation, and 
homonuclear molecules. Under these assumptions, the 
electronic bound-state wave function reads 

lfo(r) = C^ (r - R/2) + e0 o (r + R/2)), (1) 

where e = ±1, C 4 , = 1/^2(1 + eS(R), with 

S(R) = J [Mr - R/2)]* <Mr + R/2)d 3 r (2) 

The positive and negative signs for e denote symmetric 
and antisymmetric orbitals, respectively. For simplicity, 
unless otherwise stated we will consider parallel-aligned 
molecules. 

The SFA transition amplitude for high-order harmonic 
generation reads, in the specific formulation of Ref. Q 
and in atomic units, 

/>oo />t /• 

M< n > =i dt dt' / d 3 kd* ec (k(t))d ion {k(t')) 

J —OcJ — OG J 

exp[iS(t,t',Q,k)]+c.c, (3) 



£dr[k + A(r)}=0, (6) 

and 

2(0 -J p ) = [k + A(t)f. (7) 

Eq. ([5|) gives the conservation of energy at the instant 
i'of ionization, and has no real solution. Indeed, the time 
t' will possess a non- vanishing imaginary part. This is 
due to the fact that tunneling is a process which has no 
classical counterpart. In the limit I p — > 0, corresponds 
to the physical situation of a classical electron reaching 
the continuum with vanishing drift velocity. Eq. ([6]) 
expresses the fact that the electron propagates in the 
continuum from t' to t, when it returns to the site of its 
release. Eq. (0 yields the conservation of energy at the 
recombination instant t, when the kinetic energy of the 
returning electron is converted into high-order harmonic 
radiation. 

One should note that the transition amplitude ([3]) is 
gauge dependent [H, [23| • Firstly, the interaction Hamil- 
tonians H int (t'), which are present in d lon (k(t')), are dif- 
ferent in the length and velocity gauges. Furthermore, 
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in both velocity- and length-gauge formulations, field- 
free bound states are taken, which are not gauge equiv- 
alent. Therefore, different gauge choices will yield dif- 
ferent interference patterns [1, 111, H3, [13, El, Q • This 
problem has been overcome to a large extent by con- 
sidering field-dressed bound states, as a dressed state in 
the length gauge is gauge-equivalent to a field-free bound 
state in the velocity gauge, and vice- versa (for details see 



B. Double-slit interference condition 



The matrix element d roc (k) = ^k r ■ e x \tp ) is then 
given by 



4S(k) 



(2tt) 3 / 2 



■cosWdpXk) 



Rx 



sin(0)0(k) 



for bonding molecular orbitals (i.e., e > 0), or 



(8) 



4e1(k) 



(2tt)3/2 



sin(0)d Pi >(k) - ^ cos(0)^(k) 



in the antibonding case (i.e., e < 0), with $ = k-R/2. In 
the above-stated equations, R x denotes the projection of 
the internuclear distance along the direction of the laser- 
field polarization. 

In Eqs. ||5J) and ([9]), the terms with a purely trigono- 
metric dependence on the internuclear distance yield the 
double-slit condition in |5j. The maxima and minima in 
the spectra which are caused by this condition are ex- 
pected to occur for 



k • R = 2nn and k • R = (2n + l)ir, 



(10) 



respectively, for bonding molecular orbitals (i.e., e > 0). 
For antibonding orbitals, the maxima occur for the odd 
multiples of 7r and the minima for the even multiples. 
In the length and velocity gauges k(r) = k + A(r) and 
k(r) = k, where r = t, t', respectively. 

The remaining terms grow linearly with R x , and are an 
artifact of the strong-field approximation, due to the fact 
that the continuum states and the bound states are not 
orthogonal in the context of the strong-field approxima- 
tion [3, HI, Hl| • F° r that reason, they will be neglected 
here (for rigorous justifications see [1, HI]). 

In the length gauge, d TCC (k(t)) = di on (k(t')), with 
k(t) = k + A(t), while in the velocity gauge, 

«e(k) = 2 cosw(k), (id 



(2tt)3 



, (a) .-._ ^[k + Aft')] 2 



(2tt) 3 / 2 



sin(^(k), (12) 



with k(t) = k, for bonding and antibonding molecu- 
lar orbitals, respectively. The simplest and most widely 



adopted [1, H, [H, 03, d, HI] procedure is to employ 
the prefactors di on (k) and d re c(k) and the single- atom 
saddle-point equations ©-([T]). In this case, we consider 
the origin, from which the electron leaves and returns, as 
the geometric center of the molecule. 

C. Modified saddle-point equations 

The prefactors d[^ n (k) and drec(k) will now be expo- 
nentialized and incorporated in the action (for details, 
see [H, [HI)- For the recombination matrix element, we 
take the expression 



2iCV 
(2tt) 3 / 2 



f R 
cos I k • — 



(13) 



for which the spurious term is R x is absent. In the ex- 
pression for the antibonding case, the cosine term in Q13p 
should be replaced by sin(k • R/2). Without loss of gen- 
erality, the same procedure can also be applied to more 
complex orbitals. 

This leads to the sum 



3=1 f=l 



3 V 



(14) 



of the transition amplitudes 
Cij. 



M jv = 



f dt' J dt J d 3 pr](k,t 1 t') 



(2tt) 3 / 2 
x exp[iSjv(k, Q,t, t')], 



(15) 



%>(k(t)) 9g s 0(k(t')). The terms 



with r)(k,t,t') 

Sj U (k, f2, t, t') correspond to a modified action, which in- 
corporates the structure of the molecule. Explicitly, they 
read 

s ju (k, n, t, f) = s(k, n, t, t 1 ) + {-i) v+1 £(R, t, (i6) 

where £{R,t,f) = [k(t)-(-l) v+ ik(t')] • R/2. 

We will now compute the amplitudes Mj V employ- 
ing saddle-point methods. For this purpose, we will 
seek values for t, t' and k which satisfy the condi- 
tions duS jv (k,Q,t,t') = 0, d t S jv (k,Q,t,t') = and 
dt'SjuikfCl^ft') = 0. This leads to the saddle-point 
equations 



[k + A(0] 5 



-I p + {-lf v+ ' +1 d v k{t') (17) 



with Q = 



J\k + A(s)]ds + (-ir +1 d k C = 0, 



k(i) - {-iy+»k(t') ■ R/2 and 



[k + A(i)] 5 



n-I p + (-iydtk(t) • R/2. 



(18) 



(19) 
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Eq. (fT7)) corresponds to the tunnel ionization process, 
saddle-point equation (fT5|) gives the condition that the 
electron returns to its parent molecule and Eq. (fl9|) 
expresses the conservation of energy at the instant of 
recombination, in which the kinetic energy of the elec- 
tron is converted into high-order harmonic radiation. 
The above-stated saddle-point equations depend on the 
gauge, on the center Cj from which the electron was freed 
and on the center C u with which it recombines. Below 
we will have a closer look at specific cases. We will start 
by analyzing Eqs. (fTT)) and (fT9)) . which, physically, cor- 
respond to the ionization and recombination process, re- 
spectively. 

If the length gauge is chosen, both equations are ex- 
plicitly written as 



IV1 



12 



r R T 



C. " 0. 



M 



21 



R 



M 



22 



R 



[k + ^ (t,)] " = -I p + {-l) 2u+j E(t') ■ R/2, (20) 



and 



[k + A(t)] 2 



Cl-I p + (-l) v+1 E(t)-Ti/2, (21) 



FIG. 1: Schematic representation of the four possible recom- 
bination or rescattering scenarios described by Eq. (|14|l . The 
centers C\ and C2 in the molecule, as well as the transition 
amplitudes Mj v (j, v — 1, 2) are indicated in the figure. 



respectively. For this specific formulation, there exist 
potential-energy shifts on the right-hand side, which de- 
pend on the external laser field E(r)(r = t, if) and on the 
internuclear distance R. At the ionization or recombina- 
tion times, depending on the center, they increase, or sink 
the potential-energy barrier through which the electron 
must tunnel, or the energy of the state with which it will 
recombine. In the specific case discussed here, there is a 
decrease in the barrier at C2 and an increase at C\. Their 
meaning and existence altogether has raised considerable 
debate in the literature @, M, [U, Hf| . 

In the velocity gauge, the saddle-point equations (jTTJ) 
and HI]) read 



[k + A(£')] 2 



v 



(22) 



the molecule. In the former and latter case, the return 
condition reads 



J [k + A(s)]ds = 0, 



or 



J [k + A(s)]d 



s + (-1) V+1 R = 



(24) 



(25) 



In Eq. (|25p , the index v = 2 corresponds to the transition 
amplitudes M 12 (center C\ to center C 2 ) and M 2 \ (center 
C 2 to center C\), respectively. For clarity, the scenarios 
described above are summarized in Fig. 1. 



and 



[k + A(t)f 



f2 — I„ 



(23) 



These equations do not exhibit the above-mentioned 
potential-energy shifts, and resemble the saddle-point 
equations obtained for a single atom ||. Furthermore, 
if the limit I p — > is taken, Eq. (f2"2")) describes a classi- 
cal particle reaching the continuum with vanishing drift 
momentum. In contrast, in the length gauge neither the 
classical limit nor the single-atom equations are obtained. 

We will now discuss the saddle-point equation (TIE)) 
which gives the return condition. For both length and 
velocity gauges, one may distinguish two main scenar- 
ios: either the electron leaves and returns to the same 
center, i.e., v — j, or the electron is freed at a center 
Cj and recombines with the other center C v , j ^ v, in 



III. QUANTUM INTERFERENCE AND 
DIFFERENT RECOMBINATION SCENARIOS 

In the following we will discuss high-order harmonic 
spectra. For simplicity, we will consider that the elec- 
trons involved are initially bound in Is states. This gives 



<Xk) 



[k 2 



(26) 



in the high-order harmonic prefactors dj^(k) and fftec(k). 

In Fig. [21 we will commence by displaying the overall 
contributions, computed using the prefactors c?|^(k) and 

drec(k) and single-atom saddle point equations, instead 
of the modified saddle-point equations (TTT)) - (fll)j) . for the 
length and velocity gauges. For comparison, the also 
present the contribution from all transition amplitudes 
Mf?/. In the present computations, we considered up to 



5 



c 

=3 



-18 



-20 



J3 

« _ 22 

> 
o 

o 
E 

s— 
CO 

X 



-24 



-26 



-28 



-30 



I I 1 1 


i i i i i 
n= 1 






; 


_ \ 




- - Velocity gauge prefactor 




Length gauge prefactor 

i.i. 



10 20 30 40 

Harmonic Order N 



50 



60 



FIG. 2: Spectra computed employing the single-atom orbits 
and two center prefactors, for the length and velocity gauges, 
compared to the length-gauge spectrum obtained employing 
modified saddle-point equations. We consider here the mod- 
ified length form (|13[) of the dipole operator, which excludes 
the term with a linear dependence on R x . The atomic sys- 
tem was approximated by the linear combination of Is atomic 
orbitals with I p = 0.57 a.u.. The internuclear distance and 
the alignment angle are R — 2.068 a.u., and 8 = 0, respec- 
tively. The driving field intensity and frequency are given by 
/ = 3 x 10 14 W/cm 2 , and u = 0.057 a.u., respectively. The 
interference minimum at n = 1 is indicated by the vertical 
fine in the figure. The difference in the orders of magnitude 
between the velocity and length gauge spectra is due to the 
different prefactors di on (k). 



five pairs of orbits starting at the first half-cycle of the 
field, i.e., < t' < T/2. 

In the length gauge, the interference condition predicts 
interference extrema at 57 = I p + n 2 ir 2 /(2R 2 ). For the 
parameters in the figure, this yields a minimum near Q = 
31c<j, for n = 1. Even though this minimum is shallower 
if modified saddle-point equations are taken, it can be 
easily identified. 

In contrast, in the velocity gauge, the above-mentioned 
interference patterns are absent. This is due to the fact 
that the interference condition changes. The maxima and 
minima re now given by (|10|) . with k instead of k = k + 
A(t). This will lead to interference extrema at harmonic 
frequency fl = I p + [n 2 n 2 /R 2 + 2nnA(t)/R + A 2 (t)] /2. 
Roughly, if we assume that the vector potential at the 
electron return time is A(t) ~ this will correspond 

to O ~ 97lj. This frequency lies far beyond the cutoff 
( fi ~ 47cj), so that there will be a breakdown in the 
interference patterns [TJ, [26j. For this reason, in the 
following figures we will consider only the length-gauge 
situation. 



A. Modified saddle-point equations 

Subsequently, in Fig. [3l we present the contributions 
from the different recombination scenarios. In panel (a), 
the contributions from the topologically similar scenar- 
ios, involving only one or two centers, are depicted. We 
observe that the interference minimum mentioned in Fig. 
1 is absent for both types of contributions. At first 
sight, this seems to contradict the double-slit picture. In 
fact, for both \M 12 + M 2 i\ 2 and \Mn + M 22 \ 2 , high-order 
harmonic emission at spatially separated centers takes 
place. Therefore, one would expect well-defined interfer- 
ence patterns to be present. One should note, however, 
that the potential energy shifts ±E(f') • R/2 sink the po- 
tential barrier for the orbits starting at C 2 and increase 
the potential barrier for those starting at C\. Thus, the 
latter contributions are strongly suppressed and do not 
contribute significantly to the two-center interference. 

This is in agreement with panel (b), in which the con- 
tributions from the processes \Mjj + Mj V \ 2 {j, v = 1,2 
and v 7^ j) starting from the same center and ending at 
different centers are depicted. Therein, the contributions 
of the processes starting at C 2 are roughly two orders of 
magnitude larger than those starting at C\. This is due 
to the fact that the barrier through which the electron 
must tunnel in order to reach the continuum is much 
wider for the latter center. Furthermore, the two-center 
interference minimum near il = 31w is present. This is 
expected, as the contributions from the centers C\ and 
C 2 exhibit the same order of magnitude for both types 
of orbits. 

Finally, in panel (c) we display the contributions 
\Mjj + M v A 2 from the processes starting at different cen- 
ters and ending at the same center. In this case, the 
interference minimum is absent. This was expected for 
two reasons. First, for these orbits, there is no high- 
order harmonic emission taking place at spatially sepa- 
rated centers. Second, even if this were the case, the 
contributions from the orbits starting at C 2 are much 
stronger than those starting at C\ . 

Since the potential energy shifts ±E(f') • R/2 depend 
on the internuclear distance, it is legitimate to ask the 
question of whether, for small internuclear distances, a 
minimum is present in the contributions from the topo- 
logically similar scenarios. In Fig. [5J we considered such 
a situation. From the interference condition, we expect 
a minimum near f2 = 69oj. This minimum is present 
for the overall contributions, and also for the processes 



M. 



j j 



Mj U \ 2 (j,v — 1,2 and v ^ j) starting from the 



same center and ending at different centers [Fig. SJ(a)]. 
It is however absent for the interference of topologically 
similar processes [Fig. |U(b)]. This is due to fact that, 
even for this small internuclear distance, the orbits start- 
ing from C 2 lead to larger contributions than those start- 
ing from C\. Indeed, a closer look at Fig. [U(a) shows that 
the contributions \M\\ + Mi 2 \ 2 are roughly one order of 
magnitude smaller than \M 22 + M 2 i\ 2 . 

Possibly, in order to obtain well-defined maxima and 
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FIG. 3: Contributions to the high-harmonic yield from the 
quantum interference between different types of orbits, for 
internuclear distance R = 2.068 a.u. The remaining param- 
eters are the same as in Fig. 2. Panel (a): Orbits involv- 
ing similar scattering scenarios, i. e., \Mii + M22I 2 , and 
|Mi2 + M21I 2 . Panel (b): Orbits starting at the same center, 
i.e., |Mii + Mi 2 | 2 and |M 2 i+M 2 2| 2 . Panel (b): Orbits ending 
at the same center, i.e., \Mn + M 2 i\ 2 and |M i2 + M 22 | 2 . For 
comparison, the full contributions |M 2 i + M22 + AI11 +M12I 2 
are displayed as the light gray circles in the picture. The in- 
terference minimum at n = 1 is indicated as the vertical line 
in the figure. 



minima for the contributions of topologically similar sce- 
narios, it would be necessary to reduce the internuclear 
distance even more. In this case, however, none of the 
assumptions adopted in this paper, such as the LCAO 
approximation, hold. In this context, it is worth noticing 
that the parameters adopted in Fig. [4] are also somewhat 
unrealistic, as far as this specific approximation is con- 
cerned. If, however, an alternative ionization pathway is 
provided, so that the electron may reach the continuum 
without the need of overcoming the potential-energy bar- 
riers, the contributions from the topologically similar sce- 
narios may lead to well-defined patterns. Indeed, in pre- 
vious work, we employed an additional attosecond-pulse 
train in order to release the electron in the continuum, 
and obtained an interference minimum in this case 
We were, however, changing the physics of the problem 
by providing a different ionization mechanism. In the 



FIG. 4: Contributions from different types of orbits to the 
high-harmonic yield, for internuclear distance R — 1.2 a.u. 
and intensity / = 8 x 10 14 W/cm 2 . The remaining parameters 
are the same as in the previous figures. Panel (a) gives the 
contributions from the topologically similar scattering scenar- 
ios, i. e., \Mu + M22I 2 , and |Mi 2 + M 2 i| 2 , and panel (b) of 
the orbits starting at the same center, i.e., |Mn +Mi 2 | 2 and 
|M 2 i +Af 22 | 2 . 



following, we will investigate the issue of the potential- 
energy shifts for this set of parameters, employing an 
alternative method. 



B. Modified prefactors 

On the other hand, the transition amplitudes Mj v may 
also be grouped in such a way as to obtain effective pref- 
actors. Such prefactors may then be related to the quan- 
tum interference of specific types of orbits. Hence, one 
may mimic the influence of the above-stated scenarios 
even if the single-atom saddle-point equations ©-([7]) are 
taken into account. For the symmetric combination of 
atomic orbitals considered here, there would be four dif- 
ferent sets of prefactors, which are explicitly given by 

4^(k,M') = 2exp[(-l)^k(t')-R/2] (27) 
x cos[k(t) • R/2]77(k,i,i'), 
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FIG. 5: Contributions from different types of orbits to the 
high-harmonic yield, for the same parameters as in Fig. 4. 
Panel (a) gives the contributions from the topologically simi- 
lar scattering scenarios, i. e., |Mn + M22I 2 , and \Mi2+Phi\ 2 , 
and panel (b) of the orbits starting at the same center, i.e., 
[M11 + Mi 2 | 2 and A/21 + Afel 2 . All results in this figure 
have been computed mimicking the above-stated processes by 
employing the modified prefactors (27)-(30) and single-atom 
saddle-point equations. 



respectively, are involved. Interestingly, only the prefac- 
tors l ea{ i to the same interference conditions as the 
overall double-slit prefactor ((8]). 

Furthermore, one should note that, if all parameters in- 
volved were real, for the first two prefactors there would 
be the symmetry \d^ (k, t, t')\ 2 = \<$£>(k,t,lf)\* and 

|d^ 3 0M,*')| 2 = |d£ri(k,t,f)| 2 . This would lead to the 
same transition probabilities, as one transition amplitude 
is the complex conjugate of the other. This is, however, 
not the case, and can be seen by inspecting Eq. I|27[). 
Specifically in the length gauge, k(t')= k + A(t'). De- 
pending on the center, this will lead to exponentially de- 
creasing or increasing factors exp[^Im[fc + A(t')]R] in the 
transition probability \Mjj + M V j\ 2 . Clearly, this proce- 
dure is less rigorous than that adopted in the previous 
section, as we are not considering the influence of the 
potential-energy shifts in the imaginary part of t' . 

In Fig. 03 we display the results obtained following 
the above-stated procedure, for the same parameters as 
in FigHJ Once more, we see that the contributions of 
topologically similar processes, involving either one or 
two centers, do not lead to a well-defined interference 
minimum (Fig. [3(a)). Additionally, the quantum inter- 
ference of the two different kinds of processes starting 
from the same center Cj leads to a well-defined mini- 
mum at the expected frequency Q — 69to. Furthermore, 
the contributions from the orbits starting at Ci are also 
roughly one order of magnitude smaller. The main dif- 
ference between the two approaches is that the interfer- 
ence minimum is much deeper if modified prefactors are 
taken, as compared with the results obtained with mod- 
ified saddle-point equations. This discrepancy is present 
throughout, and has also been observed in Ref. jl4j . 



IV. CONCLUSIONS 



d£?(M,f) = 2exp[(-l)^[k(t)+A(t')]-R/2](28) 
x cos[k • R/2]r?(k,t,t')) 

4ame(k, t,t')=2 cos[[A(t) - A{t')} ■ R/2]7?(k, t, t% (29) 
and 

d d jff(M,0 = 2cos[p-R+[A(t)+A(i')]-R/2j;30) 
x v {k 7 t,t'). 

The prefactor corresponds to the transition am- 
plitudes Mjj + M V j in which the electron starts at the 
same center and recombines with different centers in the 
molecule. The prefactor is related to the transition 
amplitudes M V j + Mjj in which the electron starts at dif- 
ferent centers, but ends at the same center Cj. Finally, 
rfsamo and ddie corresponds to the topologically similar 
processes, in which only one, or two center scenarios, 



The results presented in this work indicate that the 
double-slit interference maxima and minima in the high- 
order harmonic spectra, which are attributed to HHG at 
spatially separated centers, are mainly due to the quan- 
tum interference between the processes \Mjj + Mj U \ 2 : 
(j = 1,2) in which the electron is released in the contin- 
uum at a center Cj in the molecule, and, subsequently, 
recombine either at the same center or at a different cen- 
ter C v . This can be seen either by employing modified 
saddle-point equations, in which the one-or two center 
scenarios are incorporated in the action, or by utilizing 
modified prefactors in which only the above-stated pro- 
cesses are included. In particular, when using the latter 
method, the transition amplitudes related to both pro- 
cesses can be grouped in such a way that the correspond- 
ing prefactor d\^ (k,t,t') exhibits the same interference 
conditions as those in the overall prefactorjf8|) . This is 
in agreement with the results obtained in [14T ]. 

These results are not obvious, as there are other pro- 
cesses which lead to high-order harmonic emission at dif- 
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ferent centers in the molecule. They do not lead, how- 
ever, to the double-slit interference patterns. This is due 
to the fact that, in the present framework, there exist 
potential-energy shifts that, depending on the center, 
sink or increase the barrier through which the electron 
must initially tunnel. Therefore, they strongly suppress 
the contributions to the spectra from one of the centers 
in the molecule. This will lead to an absence of the two- 
center interference patterns for processes starting at dif- 
ferent centers. We have verified that this suppression 
occurs even for small internuclear separations. 

Such potential-energy shifts, however, are only present 
in the length-gauge strong-field approximation and have 
raised a great deal of controversy [l(J EH EH HS|- In 
fact, it is not even clear whether they are not an arti- 



fact of the SFA. On the other hand, even if single-atom 
saddle-point equations are taken, we found a suppression 
in the yield for one of the centers of the molecule. This in 
principle counterintuitive result is related to the fact that 
the electron start time t' has a non-vanishing imaginary 
part, which suppresses or enhances the yield through the 
corresponding prefactors. 



Acknowledgments 

This work has been financed by the UK EPSRC (Ad- 
vanced Fellowship, Grant no. EP/D07309X/1). 



[1] See, e.g., A. Scrinzi, M. Y. Ivanov, R. Kienberger, and 

D. M. Villeneuve, J. Phys. B 39, Rl (2006) for a review 

onn the subject. 
[2] P. B. Corkum, Phys. Rev. Lett. 71, 1994 (1993); K. C. 

Kulander, K. J. Schafer, and J. L. Krause in: B. Piraux et 

al. eds., Proceedings of the SILAP conference, (Plenum, 

New York, 1993). 
[3] M. Lewenstein, Ph. Balcou, M. Yu. Ivanov, A. L'Huillier 

and P. B. Corkum, Phys. Rev. A 49, 2117 (1994). 
[4] W. Becker, S. Long, and J. K. Mclver, Phys. Rev. A 41, 

4112 (1990); ibid. 50, 1540 (1994); M. Lewenstein, K. C. 

Kulander, K. J. Schafer and Ph. Bucksbaum, Phys. Rev. 

A 51, 1495 (1995). 
[5] M. Lein, N. Hay, R. Velotta, J. P. Marangos, and P. L. 

Knight, Phys. Rev. Lett. 88, 183903 (2002); Phys. Rev. 

A 66, 023805 (2002); M. Spanner, O. Smirnova, P. B. 

Corkum and M. Y. Ivanov, J. Phys. B 37, L243 (2004). 
[6] G. Lagmago Kamta and A. D. Bandrauk, Phys. Rev. A 

70, 011404 (2004); ibid. 71, 053407 (2005). 
[7] V. I. Usachenko, P. E. Pyak, and Shih-I Chu, Laser Phys. 

16, 1326 (2006). 
[8] H. Hetzheim, C. Figueira de Morisson Faria, and W. 

Becker Phys. Rev. A 76, 023418 (2007). 
[9] D. B. Milosevic, Phys. Rev. A 74, 063404 (2006). 
[10] W. Becker, J. Chen, S.G. Chen, and D. B. Milosevic, 

Phys Rev. A 76, 033403 (2007). 
[11] M. Busuladzic, A. Gazibegovic-Busuladzic, D. B. 

Milosevic, and W. Becker, Phys. Rev. A 78, 033412 

(2008); Phys. Rev. Lett. 100, 203003 (2008). 
[12] R. Kopold, W. Becker and M. Kleber, Phys. Rev. A 58, 

4022 (1998). 

[13] C. C. Chirila and M. Lein, Phys. Rev. A 73, 023410 
(2006). 

[14] C. Figueira de Morisson Faria, Phys. Rev. A 76, 043407 



(2007). 

[15] C. Figueira de Morisson Faria, arXiv: 0807.2763 [atom- 
ph] 

[16] J. Muth-Bohm, A. Becker, and F. H. M. Faisal, Phys. 

Rev. Lett. 85, 2280 (2000); A. Jaron-Becker, A. Becker, 

and F. H. M. Faisal, Phys. Rev. A 69, 023410 (2004); A. 

Requate, A. Becker and F. H. M. Faisal, Phys. Rev. A 

73, 033406 (2006). 
[17] T. K. Kjeldsen and L. B. Madsen, Phys. Rev. A 71, 

023411 (2005); Phys. Rev. Lett. 95, 073004 (2005); C. 

P. J. Martiny and L. B. Madsen, Phys. Rev. Lett. 97, 

093001 (2006); C. B. Madsen and L. B. Madsen, Phys. 

Rev. A 74, 023403 (2006). 
[18] V. I. Usachenko, and S. I. Chu, Phys. Rev. A 71, 063410 

(2005). 

[19] X. Zhou, X. M. Tong, Z. X. Zhao and C. D. Lin, Phys. 

Rev. A 71, 061801(R) (2005); ibid. 72, 033412 (2005). 
[20] C. C. Chirila and M. Lein, J. Mod. Opt. 54, 1039 (2007). 
[21] C. Figueira de Morisson Faria, T. Shaaran, X. Liu and 

W. Yang, Phys. Rev. A 78, 043407 (2008). 
[22] P. Salieres, B. Carre, L. LeDeroff, F. Grasbon, G. 

G. Paulus, H. Walther, R. Kopold, W. Becker, D. B. 

Milosevic, A. Sanpera and M. Lewenstein, Science 292, 

902 (2001). 

[23] A. Fring, V. Kostrykin and R. Schrader, J. Phys. B. 29, 
5651 (1996); D. Bauer, D. B. Milosevic and W. Becker, 
Phys. Rev. A 72, 023415 (2005). 

[24] O. Smirnova, M. Spanner and M. Ivanov, J. Phys. B 39, 
S307 (2006). 

[25] C. Figueira de Morisson Faria, H. Schomerus and W. 

Becker, Phys. Rev. A 66, 043413 (2002). 
[26] O. Smirnova, M. Spanner and M. Ivanov, J. Mod. Opt. 

54, 1019 (2007). 



